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Unveiling the nuclear structure at small x using azimuthal 
correlations in dA collisions 

Kirill TuCHiN 

Department of Physics and Astronomy, Iowa State University, Ames, lA 50011 

We discuss azimuthal correlations in double inclusive hadron production in high energy 
dA collisions at RHIC. We argue that the leading logarithmic approximation is inadequate 
for description of experimental data. Realistic shape of the azimuthal correlation function is 
obtain only if we keep terms that are finite in the rapidity difference between the produced 
hadrons. Likewise, transverse momentum dependence of parton distribution functions in 
both deuteron and nucleus must be retained. We observe depletion of the back-to-back 
T-H correlations in central dAu collisions in the forward direction consistent with the Color 

Glass Condensate predictions. 

(N 

O §1. Introduction 

(D 

Q Azimuthal correlations is an important tool to investigate properties of the new 

0^ form of nuclear matter created in high energy dA and AA collisions at RHIC - 

Color Glass Condensate (CGC) - a coherent quasi-classical state of fast gluons and 
I— I quarks. In this article we report on a recent study of azimuthal correlations in dAu 

collisionJil and demonstrate how the CGC modifies the correlation function. 

In a pioneering papeJD it was proposed to study azimuthal correlations of 
hadrons produced at large rapidity separation Ay S> 1. The idea is that such 
correlations are mediated by the BFKL Pomeron - a collective state of a large num- 
ber of gluons, which is a most abandon excitation in the CGC. Therefore, unlike the 
1—1 hadron production in hard collisions, where there is a strong back-to-back correlation 

^ at opening azimuthal angle Acp = vr, correlations in the CGC should be significantly 

reduced. Although only qualitative arguments were presented iiiH' it seemed that 
^SJ a more quantitative analyses is not too difficult since Ay ^ 1 corresponds to the 

well-studied multi-regge-kinematics (MRK) regime and the production amplitude is 
CN essentially the real part of the leading order BFKL. We will show however, that the 

MRK is not sufficient to describe the data at Ay = 3 at RHIC and terms finite in 
Ay must be kept. 

It has been suggested irll' that correlations at small Ay in the forward direction 
can also be used to study CGC. Indeed, forward direction correspond to small x of 
rS nucleus where the CGC effects are strongest. They reduce both single and double 

inclusive hadron production and thus back-to-back correlations are suppressed. The 
problem is that the MRK approximation is not applicable at all in this case. One 
therefore has to rely on phenomenological models, which offer descriptions that are 
analytically accurate only in parts of the relevant kinematic region. There are two 
such approaches: one that is based on the dipole mo and another one that is 

based on the feT^-factorization.^' 

In an approrach based on the 'dipole model' one reduces the relevant scattering 
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amplitudes to a product of light-cone "wave-functions" and combinations of CGC 
field correlators in the configuration space. These correlators satisfy a set of evo- 
lution equations with certain initial conditions. The advantage of this approach is 
that it rests on accurate theoretical treatment of the gluon saturation region. In 
this approach the double inclusive gluonlD quark-anti-quarlJlHll' and valence quark 
- gluorl^f cross sections were calculated. In most cases in order to obtain analytical 
results an assumption of large rapidity gap Ay ^ 1 must be made. Only in the very 
forward region dominated by valence quarks of deuteron, i.e. at Xp ~ 1 finite Ay 
effects have been taken into account^''A' transverse momentum dependence of the 
valence quark distribution is neglected in these papers (resulting in collinear factor- 
ization of the deuteron parton distributions). Transverse momentum dependence of 
valence quark distribution seems to be very important for the shape of the correla- 
tion function, though not as important for generation of suppression which is mostly 
due to large momentum flow from the nucleus. 

Another approach is '/c-r-factorization', which assumes that 2 — )• n process and 
the two-point correlation functions of CGC fields can be factored out. In this ap- 
proximation, the 2 — )• 4 amplitudes were calculated for an arbitrary Ay (quasi multi- 
Regge kinematics, QMRK) iiJlHI3'IIll for gg — )- ggqq and iillll'lll''lll' for gg — )- gggg 
processes. Although generally fcr-factorization fails in the gluon saturation region, 
there are valid reasons to believe that it provides a reasonable approximation of the 
observed quantities. Indeed, it was proved that fcr-factorization provides the exact 
result for the cross section for s ingl e inclusive gluon production in the leading log- 



arithmic approximation (LLA) (2-4)^ (though there is a subtlety in the definition 
of the unintegrated gluon distribution (/ilSI'llIf ). Although A^T-factorization fails for 
the double-inclusive heavy quark production, the deviation from the exact results is 
not large at RHIC energiesU^f At transverse momenta of produced particles much 
larger than Qg, /cT-factorization rapidly converges to the exact results. There are 
also numerous indications that /c^-factorization is phenomenologically reliable (sed3 
for examples). 

§2. Correlations at \yT — VaI ^ 1 

First, we would like to consider correlations at small rapidity separations. Az- 
imuthal correlation function is defined as 

C{AcP) = J-J^ (2-1) 

where dN/d{A(j)) is the number of pairs produced in the given opening angle Acj) 
and A'^trig is the number of trigger particles. The number of pairs is given by 

dN ./■.../■. /",../■./ dA-trig dA^ass dA„ 



27r / dhrkr / dyT / dkA^A / dyA ,o, , ,n, , h 



d{A(j)) J J J J yd'^kxdyT d'^kAdyA d^kxdyx d^kAdyA 

(2-2) 

where k^ and yx are the transverse momentum and rapidity of the trigger particle 
and k^ and yA are the transverse momentum and rapidity of the associate one. 
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We denote kx = etc. throughout this paper. The first term on the r.h.s. of 



(2-2) corresponds to gluon production in two difi^erent sub-cohisions (i.e. at difi'erent 
impact parameters) and therefore gives a constant contribution to the correlation 
function, whereas the second term on the r.h.s. describes production of two particles 
in the same sub-collision. The number of the trigger particles is given by 

N,.,, = 2,JikrkrJdyr^^. (2-3) 

Expression for the single inclusive gluon cross section is well-known (see e.gEZl). The 
corresponding multiplicity reads 



In the center-of-mass frame x-t = -7= explztyj . Equation (2-4) is derived in the 
multi-Regge kinematics (MRK) x± -^1. 

The correlated part of double-inclusive parton multiplicity is given by 



f d^qi f d'q2 r2( ^ 11 



d^kx dyx d'^kAdyA it'^CfS± J q\ J 

X ipD{xi,ql) ipA{x2,ql) A{(i^,(i^,\T^'^A^yT - Va) , (2-5) 

where xi^2 = (fc^e^^^ + kAe^^'^)/y/s. The amplitude A was computed in the quasi- 
multi-Regge-kinematics (QMRK) iiJllUllf and recently re-derived i rJlD (the gg 
ggqq part was calculated before iiJU'^Iil'). In QMRK one assumes that xi,X2 ^ 1, 
but Ay is finite. Explicit expression for A can be found inllll 

For numerical calculations we need a model for the unintegrated gluon distribu- 
tion function (p. In spirit of the KLN mo de© we write 

^{x,q') = ^,^{1- e-«^/«^) (1 - x^ . (2-6) 
zvr^ as 

where the saturation scale of nucleus is = A^^^Q^p, with Q^^ the saturation scale 
of proton fixed by fits of the DIS data. The coupling constant is fixed at as = 0.3. 

It has been pointed out inlli' that due to 1 — 2 gluon splittings the double- 
inclusive cross section has a collinear singularity at s — t- 0, i.e. it is proportional 
to [{Ay)"^ + (Z\0)^]^^. Such singularities are usually cured at the higher orders of 
the perturbation theory. Additional contributions to the small angle correlations 
arise from various soft processes including resonance decays, hadronization, HBT 
correlations etc. Because the small angle correlations are beyond the focus of the 
present paper we simply regulate it by imposing a cutoff on the minimal possible 
value of the invariant mass s. This is done by redefining the amplitude as ^ — t- 
As/{n'^ + s). For each kinematic region, parameter ^ is fixed in such a way as to 
reproduce the value of the correlation function in pp collisions at zero opening angle 
Acj) = 0. 
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/cT^-factorization is known to give results that are in qualitative agreement with 
a more accurate approaches, but misses the overall normalization. Therefore, in 
order to correct the overall normalization of the cross sections we multiply the single 
inclusive cross section (|2-4| by a constant Ki and the double- inclusive one (2-5) by a 
different constant K22^'^ The correlation function C depends on both Ki and K2. 
However, the difference Ca = C{A(p)—C{A(l)o) depends only on the ratio K2/K1. We 
choose A(j)o in such a way that C{A(f)Q) is the minimum of the correlation function. 
This is analogous to the experimental procedure of removing the pedestalMI The 
overall normalization of the correlation function K2 / Ki - which is the only essential 
free parameter of our model - is fixed to reproduce the height of the correlation 
function in pp collisions. 

The results of the numerical calculations are shown in Fig. [l}Fig. [2] and Fig. [Sj 
In these figures we observe suppression of the bak-to-back correlation in dAu as 
compared to the pp ones in agreement with the experimental data. In Fig. [3] we also 
see the depletion of the back-to-back correlation as a function of centrality. Note, 
that at the time of publication the precise centrality classes of the data shown in the 
lower row of Fig. |3]were not known. 




Fig. 1. Correlation function at the central rapidity. Kinematic region is 4 < pT < 6, 2 < < Pt 
(all momenta are in GeV), j/t = 3.1, j/a = 3. Left (right) panel: minbias pp {dAu) collisions. 
Data fromHH 



In addition to gg — gggg and gg — ggqq processes that we took into account 
in this section, production of valence quark of deuteron gq^ — )• gqvgg gives a sizable 
contribution at forward rapidities due to not very small value of x associated with 
deuteron (x ^ 0.2 for px = 2 GeV at 7/ = 3). Contribution of this process to az- 
imuthal correlations was analyzed iiHl' in the framework of the dipole model in MRK. 
However, the corresponding expression in fc^-factorization in QMRK is presently un- 
known thus preventing us from taking it into account in our calculation. In-spite of 
this we believe that the general structure of the correlation function as well as its 
centrality dependence is not strongly affected by the valence quark contribution. We 
plan to address this problem elsewhere. 

So far we neglected the effect of fragmentation. Convoluting the correspond- 
ing cross sections with the fragmentation functions taken froniHll' results in minor 
modifications of azimuthal distributions. For example, modification of the correla- 
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Fig. 3. Correlation function at forward rapidities. Kinematic region is pr > 2, 1.5 < pa < Pt (all 
momenta are in GeV), yT = 3.1, yA = 3. Upper left (right) panel: minbias pp [dAu) collisions. 
Lower left (right) panel: peripheral (central) dAu collisions. Note: centrality of the theoretical 
calculation may not coincide with the centrality of the data (the former is not yet known at the 
time of publication). Data fromHH 



tion function of Fig. [T] is exhibited in Fig. |4j Similar conclusion holds also for the 
correlation functions at forward rapidity. 
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Fig. 4. Effect of fragmentation on the azimuthal correlation function. Solid lines are the same as in 
Fig. [l] Dashed lines represent a conservative estimate of the fragmentation effect as discussed 
in the text. 



§3. Correlations at \yT — UaI ^ 1 

If the trigger and associated particles are well separated in rapidity so that 
Iut — VaI ~ l/os ^ 1, we can apply the MRK approximation {Ay — )• oo) to the 
double-inclusive cross section. The result then factorizes into a product of two ladder 
rungs each given by the real part of the LO BFKL kernel. The corresponding formula 
is 



dNcarv Nr.al 1 



j (fqi LpD{xi,ql)ipA{x2,{'^T + '^A- 



cPkT duT (PkAdyA tt'^CpS± k^k\ 

(3-1) 

The advantage of the MRK approximation is that it allows taking into account a 
possible multi-gluon production in the interval between yx and y^- Unfortunately, in 
this approximation one also looses many features of the azimuthal angle dependence 
that are important for description of the backward pick at intermediate Ay, see 
Sec. [2j Thus, we are facing a dilemma: either to use formulas of Sec. [2] that give 
precise dependence on Ay but neglect evolution in the gap, or to take into account 
the evolution as discussed below but in the MRK limit Ay — )• oo. At present, there is 
no approach that would interpolate between these limits at intermediate Ay relevant 
for RHIC. Therefore, in this section we will calculate the correlation function in two 
limits, compare our results with the data and try to learn which approximation is 
more phenomenologically important at Ay = 3. 



Eq. (3-1) does not take into account a possible gluon emission in the rapidity 
interval between yr and yA- This is important when \yi — y2\ > '^/oig and may be 
important for the experimentally measured forward-backward rapidity correlations. 
Evolution in between the rapidities of the produced particles can be included using 
the the AGK cutting rules and the known properties of the BFKL equation as 



d'^kx dyx d'^kAdyA tt'^CfSx k'^n^ 



X G(|q^ - ky|, Iq^ - k^|,yT - y^) , (3-2) 
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where G is the Green's function of the BFKL equation. It can be written as 



G{qi,q2,y) = ^2cos(nqi • ci2)Gn{qi,q2,y) 



(3-3) 



n=0 



where functions G„ are can be found elsewhere (see e.gEl'). 

It can be easily seen that G « 2[Go+Gi cos(qi-q2)] is a verygood approximation 
at Ay = 30' Numerical calculations are presented in Fig. [s] together with the 
experimental data from UDjM I We observe that the shape of the correlation function 
is better described by (2-5) in agreement with the observation oil^'llil' that finite 



Ay corrections to the MRK approximation are essential for description of azimuthal 
correlations. On the other hand, it seems that to explain the magnitude of depletion 
one also needs to include the small x evolution effects in the gap Ay = yr — VA- 
Obviously, a more accurate description requires additional theoretical investigation of 
the finite Ay corrections. Data fronl^''^ also shows a significant isospin effect (not 
displayed here) that probably originates in the valence quark contribution not taken 
into account in the present work. This isospin effect obscures the CGC contribution 
and requires a detailed analyses that we plan to do elsewhere. At LHC one can get 
rid of the isospin effect by considering correlations at large rapidity gaps away from 
the fragmentation regions. 




Fig. 5. Forward-backward correlations. Kinematic region is pr > 2, 1 < pA < Pt (all momenta 
are in GeV), yr — 3, i/a ~ 0. Left (right) panel: minbias pp (dAu) collisions. Solid lines: 



calculations with (2-5 1 (exact 2—^4 amplitude, no evolution between the trigger and associate 



particles). Dashed line: calculations with ( |3-2[ ) (MRK approximation of 2 — >■ 4 amplitude, 
includes evolution between the trigger and associate particles). Data frorrElUlll (forward 7r° 
and midrapidity h^). 



§4. Conclusions 

We calculated the azimuthal correlation function in dAu collisions using the 
approach developed by us beforeIl2f'll3 The results are presented in Figs, 
and [5} We demonstrated that CGC is responsible for depletion of the back-to- 
back correlations in dAu collisions as compared to those in pp ones at small rapidity 
separations - at midrapidity and forward rapidity - and at large rapidity separations. 
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Our results quantitatively confirm earlier arguments ofH) 
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